AD-A11J  542  COLO  REGIONS  RESEARCH  AND  ENGINEERING  LAB  HANOVER  NH  F/G  8/12 

SEA  ICE  DRAG  LAWS  AND  SIMPLE  BOUNDARY  LAYER  CONCEPTS.  INCLUDING— ETC <U> 
FEB  82  M  G  MCPHEE 

UNCLASSIFIED  CRREL-82-4  nL 


« 


uu-ULSLt, 

REPORT  82-4 

^  S‘mp,e 

application  to  nw%%™ludlnS 


site?* 

Crtlsi  _ ~ 


O) 

tO 

00 


For  conversion  of  SI  metric  units  to  U.S./ 
British  customary  units  of  measurement 
consult  ASTM  Standard  E380,  Metric  Prac¬ 
tice  Guide,  published  by  the  American  Socie¬ 
ty  for  Testing  and  Materials,  1916  Race  St., 
Philadelphia,  Pa.  19103. 


Cover:  Melting  ice  floes  in  the  Weddell  Sea 
after  advection  into  above-freezing 
ocean  surface  water. 


CRREL  Report  82-4 

February  1982 


Sea  ice  drag  laws  and  simple 
boundary  layer  concepts,  including 
application  to  rapid  melting 


Miles  G.  McPhee 


Prepared  I  Of 

OFFICE  OF  NAVAL  RESEARCH 


Approved  tar  puttie  releeae:  dtatrlbutlon  unlimited. 


SECURITY  CLASSIFICATION  OP  THIS  PAGE  (Wlitl  Onto  Bntond) 


REPORT  DOCUMENTATION  PAGE 


iwt  jjvurTTrri 


CRREL  Report  82-4 


«.  TITLE  (mid  Subtitle) 

SEA  ICE  DRAG  LAWS  AND  SIMPLE  BOUNDARY  LAYER  CON¬ 
CEPTS,  INCLUDING  APPLICATION  TO  RAPID  MELTING 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CATALOG  NUMBER 

/Z- 


S.  TYRE  OP  REPORT  A  PERIOD  COVERED 


S.  PERFORMING  ORG.  REPORT  NUMBER 


7.  AUTHOR^) 

Miles  G.  McPhee 


e.  contract  or  grant  number^ 
Office  of  Naval  Research, 

MIPR  N0OOM80MPOOOO4 


».  PERPORMING  ORGANIZATION  NAME  ANO  ADDRESS 

U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
Hanover,  New  Hampshire  037SS 


II.  CONTROLLING  OPPICE  NAME  ANO  AODRESS 

Office  of  Naval  Research 
Arlington,  Virginia  22217 


10.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

February  1982 


IS.  NUMBER  OP  PAGES 

25 


i 


MONITORING  AGENCY  NAME  a  ADORESSfif  dffforanl  I 


i  Controlling  Ollleo)  I  IS.  SECURITY  CLASS,  (ot  tttlm  report) 


Unclassified 


Bn.  OECLASSIPICATION/DOWNGRADING 

schedule 


IA.  DISTRIBUTION  STATEMENT  (o I  Nil*  Roport) 


Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (ot  Ike  nbotrnet  ontorog  In  Block  20,  It  different  from  Report) 


19.  KEY  WORDS  (Continue  on  too 

Boundary  layer  flow 

Drag 

Drift  / 


woo  oldo  II  nocoooory  end  Identity  by  block 

Melting. 

Pack  ice 
Sea  ice 


_  BSTRACT  (ComUnuo  more  wore*  otdm  H  nowoeooty  omg  Identity  k y  block  number) 

Several  proposed  methods  for  treating  the  momentum  flux  between  drifting  sea  ice  and  the  underlying  ocean  are 
interpreted  in  terms  of  simple  planetary -boundary -layer  (PBL)  turbulence  theory.  The  classical  two-layer  approach,  in 
which  the  solution  for  a  thin  surface  layer  is  matched  to  an  Ekman  solution  for  the  outer  layer,  is  used  to  derive  several 
forms  for  the  drag  law.  These  forms  range  from  linear  (where  stress  is  proportional  to  relative  speed),  through  quadratic 
(whare  strtsa  ia  proportional  to  relative  speed  squared),  to  a  Rossby-similarity  law  like  that  used  to  express  frictional 
drag  on  geoitrophic  wind  in  the  atmosphere.  Only  formulations  which  conform  with  Rossby-similarity  scaling  are 
conaistent  with  free-drift  data  from  the  1975  AIDJEX  drift  station  experiment.  We  show  how  a  two-layer  model,  id 
which  the  non  dimensional  stress  is  allowed  to  vary  realistically  through  a  surfesAgycr  of  constant  non  dimensional 


oo  .jSTn  un  EornoH  or  i  mov  as  is 


SECURITY  CLASSIFICATION  OP  THIS  RACE  | 


i  Data  Mntorod) 


mcuwiTv  ctwtriCATiow  or  thi»  ww*; Data  inwio 

thickness,  provides  an  analytic  solution  for  the  steady-state  PBL  equation  quite  similar  to  recent  numerical  solutions. 
The  theory  is  extended  to  include  drag  reduction  due  to  buoyancy  from  rapid  melting  and  is  shown  to  agree  with 
atmospheric  results  for  geostrophic  drag  under  analogous  conditions  of  radiational  cooling.  The  theory  provides  a  basis 
for  estimating  trajectories  and  melt  rates  of  floes  drifting  into  water  warmer  than  the  ice  melting  temperature. 


Accession  "or  > 

"ifis  9* '.AT  0 

otic  n  □ 

Just if lcatlcn- - 

*T— - - - — 

Dlatrlbut Ion/ 

Availability  Codon 
jXvail  and 'or 
Dlst  Special 


S1CURITV  CLASSIFICATION  OF  THIS  PAOtfWtlF"  Data  KMatad) 


PREFACE 
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Analysis  of  the  Physical  Oceanographic  Data  Collected  on  FRAM  /. 
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viewing  the  manuscript  of  this  report. 


CONTENTS 


Abstract.. 


Background . 

Hierarchy  of  drag  laws  and  simple  models . 

The  momentum  equation  for  the  planetary  boundary  layer.. 

linear  eddy  viscosity  —  the  constant  stress  layer . 

Two-layer  eddy  viscosity . 

PBL  scaling . 

A  dimensionless  two-layer  system . 

A  dimensionless  two-layer  system  with  modified  stress . 

Evaluating  the  drag  laws . 

Rossby  similarity  parameters  and  buoyancy  effects . 

Discussion . . 

Literature  cited . 


ILLUSTRATIONS 


Kinematic  stress  vs  relative  speed  between  ice  and  ocean  under  free-drift 

conditions . 

Schematic  of  friction  velocity/surface  velocity  relation  in  complex  plane.. 

Diagram  of  turbulent  flux  mechanism . 

Two-layer  system  with  constant  eddy  viscosity . 

Two-layer  system  with  linear  eddy  viscosity  near  surface . 

Two-layer  system  nondimensionalized  by  Rossby-similarity  scales . 

Two-layer  system  in  which  stress  is  allowed  to  vary  through  a  logarithmic 

surface  layer . 

Nondimensional  stress  and  velocity  in  the  planetary  boundary  layer . 

Comparison  of  six  drag  law  formulations . 

The  Rowby -similarity  parameters  >4  and  5  as  functions  ofu,  =  ML, . 

Data  from  Clarke  and  Hess . 


TABLE 


Table 

1 .  Summary  of  drag  laws . 


SEA  ICE  DRAG  LAWS  AND  SIMPLE 
BOUNDARY  LAYER  CONCEPTS, 
INCLUDING  APPLICATION  TO 
RAPID  MELTING 


Mites  G.  McPhee 


BACKGROUND 

The  drift  of  sea  ice  is  determined  by  a  balance 
of  forces  including  a  flux  of  momentum  from  the 
ice  to  the  ocean,  which  is  more  often  than  not  of 
the  same  order  of  magnitude  as  the  tangential  stress 
exerted  on  the  ice  by  the  wind.  Its  representation 
is  thus  a  major  element  in  a  mathematical  formula¬ 
tion  of  ice-drift  dynamics,  and  it  has  received  much 
attention  as  a  problem  in  planetary-boundary-layer 
(PBL)  theory,  beginning  with  Ekman  (1905).  Most 
ice  modelers  are  reluctant  to  include  complete  sub¬ 
models  of  the  oceanic  boundary  layer  in  their  nu¬ 
merical  codes  and  seek  instead  to  parameterize  the 
tangential  stress  between  the  ice  and  ocean  as  a 
function  of  the  ice  velocity  relative  to  the  undis¬ 
turbed  ocean  by  using  relatively  simple  expressions 
which  we  shall  call  oceanic  drag  laws.  Several  drag 
laws  have  been  proposed  by  various  authors,  and  a 
modeler  perhaps  not  well-versed  in  PBL  theory  faces 
the  task  of  choosing  a  formulation  suited  to  his 
work.  The  primary  intent  of  this  report  is  to  assist 
in  that  choice  by  providing  a  framework  in  which 
several  drag  laws  may  be  compared,  both  in  terms 
of  how  well  they  describe  pertinent  data,  and  how 
they  fit  within  the  guidelines  of  PBL  theory. 

Before  the  PBL  equations  could  be  routinely 
solved  numerically,  drag  laws  for  sea  ice  (and  for  the 
atmosphere)  were  obtained  by  patching  an  Ekman- 
IJce  solution,  appropriate  for  the  outer  part  of  the 
boundary  layer,  to  the  solution  for  a  relatively  thin 
layer  adjacent  to  the  surface  (e.g.  Rossby  and  Mont¬ 
gomery  193S).  We  shall  adopt  the  same  approach, 
at  the  same  time  attempting  to  rationalize  it  with 
recent  measurements  by  progressively  introducing 
concepts  from  modem  turbulence  theory  into  a 
series  of  boundary  layer  formulations  which  provide 
drag  laws  ranging  from  tinear  (u.  where  stress  is 
proportional  to  relative  speed),  through  quadratic 
(where  stress  is  proportional  to  relative  speed  squared), 
arriving  Anally  at  a  Rossby-similarity  law  consistent 


with  one  we  proposed  earlier  on  the  basis  of  numer¬ 
ical  modeling  (McPhee  1979).  Most  drag  laws  found 
in  the  literature  (Shuleikin  1938,  Reed  and  Camp¬ 
bell  1962,  Rothrock  1975,  McPhee  1975  and  1979, 
etc.)  fall  into  one  of  the  categories  treated. 

In  addition,  the  theory  is  extended  to  include 
the  effects  of  rapid  melting  that  might  occur  if  sea 
ice  were  advected  into  comparatively  warm  water. 

The  problem  is  addressed  both  for  the  practical 
purpose  of  interpreting  and  predicting  sea  ice  be¬ 
havior  near  the  ice  margin  and  also  as  a  novel  appli¬ 
cation  of  stable  atmospheric  boundary  layer  studies. 

Despite  the  simplicity  of  the  concepts  and  math¬ 
ematical  models  developed  here  (nearly  all  the  cal¬ 
culations  were  done  with  a  hand-held  calculator- 
albeit  one  designed  for  easy  manipulation  of  complex 
numbers),  the  final  result  is  remarkably  consistent 
with  the  knowledge  we  have  of  momentum  transfer 
in  a  neutral  or  stably  stratified  PBL.  Thus  a  second 
purpose  of  this  report  is  to  emphasize  that  the  ice 
drift  problem  can  be  turned  around  to  reveal  much 
about  how  boundary  layers  work.  We  took  this 
approach,  for  example,  in  using  ice  drift,  surface 
stress,  and  current  meter  measurements  to  corrobor¬ 
ate  results  from  a  numerical  PBL  model  (McPhee 
1979),  and  also  to  investigate  time-dependent  re¬ 
sponse  to  rapidly  changing  winds  (McPhee  1980a). 

Here  we  bypass  some  of  the  more  complicated  features 
explored  numerically  in  the  earlier  works  (e.g.  inertial 
waves  and  abrupt  changes  in  density  at  the  base  of 
the  mixed  layer)  to  examine  what  the  stress/velocity 
relationship  for  drifting  sea  ice  implies  about  the 
basic  structure  of  PBL  turbulence. 

To  start  we  establish  some  conventions,  the  first 
being  that  horizontal  (two-dimensional)  vectors  are 
represented  as  complex  numbers  such  that  A  B  Axex 
+ Ayty  (where  ex  and  ey  are  orthogonal  unit  vectors) 
becomes  A  *  A  (cos  a  +  f  sin  a)  *  Ae* ,  where  i1  *  -1 . 
The  real  axis  is  aligned  with  ex,  the  imaginary  axis 
with  ex,a*  tan"1  (Ay/Ax)  and  A*  Ul .  Note  that 
a  symbol  with  a  caret  indicates  a  complex  (vector) 


quantity,  while  the  same  symbol  without  a  caret 
denotes  only  magnitude.  With  this  convention,  it 
is  meaningful  to  describe  products  or  quotients  of 
vectors,  which  are  also  vectors,  e.g.  for  B  *  Be#, 

A/B  =  (AIB)e«*-<». 

Note  also  that  rotation  of  a  vector  in  the  horizontal 
plane  through  an  angle  0  (positive  counterclock¬ 
wise)  occurs  if  the  vector  is  multiplied  by  the  factor 
eie .  Complex  notation  is  introduced  both  to  avoid 
intricate  geometrical  constructions  and  as  a  compu¬ 
tational  aid. 

The  second  convention  is  that,  unless  otherwise 
noted,  the  velocity  appropriate  for  a  specific  drag 
law  is  the  relative  velocity  between  the  ice  and  the 
undisturbed  ocean.  This  relative  velocity  is  the  ice 
velocity  with  respect  to  water  below  the  maximum 
penetration  of  the  frictional  PBL,  which  is  also  the 
geostrophic  flow  due  to  sea  surface  tilt  (provided 
there  are  no  horizontal  density  gradients  in  the 
mixed  layer).  Usually,  the  ice  drift  is  larger  than 
geostrophic  flow  and  the  momentum  flux  is  down¬ 
ward,  but  the  flux  direction  is  unimportant  as  long 
as  the  drag  law  is  expressed  in  terms  of  relative 
velocity. 

With  complex  notation  in  a  frame  of  reference 
fixed  to  earth,  the  balance  of  forces  acting  at  a 
point  in  a  large,  uniform  slab  of  ice  in  steady  drift  is 

=  p,r,  -pwrQ  +  J  +  S.  (1) 

where  m  -  the  ice  mass  per  unit  surface  area 
/=  the  Coriolis  parameter 
Vt  -  the  absolute  ice  velocity 
p,  and  pw  =  air  and  water  densities 
t,  3  the  kinematic  wind  stress 
ro  3  the  kinematic  stress  exerted  by  the 
ice  on  the  ocean 

/*  the  force  due  to  a  gradient  of  internal 
ice  stress 

S  3  a  gravitational  force  due  to  sea-surface 
tilt  (S  3  -mgtf) 

f  3  the  elevation  of  the  sea  surface 
g  3  the  acceleration  of  gravity. 

.  By  definition,  geostrophic  flow  in  the  ocean 
Vt  is  given  by 

-got*  Sim  (2) 

so  that 

mfiyt  -  Ft>e"/*  3  imfV=  p,T,  -  pwrG  ♦  / 
where  V  i*  the  relative  velocity  described  above. 


In  McPhee  (1979)  we  analyzed  the  relationship 
between  tw  and  V,  using  wind  stress  and  drift  velocity 
measured  during  the  197S  summer  melt  season  at 
four  drifting  stations  of  the  AtDJEX  Experiment 
near  the  center  of  the  Beaufort  Gyre  in  the  centra] 
Arctic.  The  major  difference  between  this  and 
earlier  analyses  (the  most  extensive  being  that  of 
Reed  and  Campbell  1962)  is  that  we  considered 
only  data  from  a  60-day  period  during  which  we 
were  reasonably  confident  that  the  internal  ice  force 
/  was  small.  During  the  remainder  of  the  year-long 
AIDJEX  drift,  /  was  often  a  significant  part  of  the 
force  balance  (McPhee  1980b),  and  affected  the 
long-term  statistics  of  ice  drift  appreciably.  During 
the  short  time  of  the  year  when  /  is  negligible, 
measurements  of  V%  and  the  1 0-m  wind  velocity, 

Ui0,  along  with  estimates  of  Kg,  m,  and  the  wind 
drag  coefficient,  C)  0,  provide  the  interfacial  stress 
to  from  the  relation 

Pw fD  =  p,C, oUioUio  - itnfV ■  (3) 

Given  a  sufficient  population  of  sample  pairs,  rD 
and  V,  an  expression  of  the  form 

r0=aV* 

can  be  derived  statistically  by  a  linear  regression  of 
the  equation 

IogT0  =  4  log  V  +  log  a. 

Figure  1  from  McPhee  (1979)  shows  results  of  such 
an  analysis  on  summer  AIDJEX  data  filtered  to 


ICE  SPEED  (cm  »'') 

Figure  1.  Kinematic  stress  vs  relative  speed  between 
ice  and  ocean  under  free-drift  conditions.  Solid  curve 
h  best-fit  exponential,  dashed  curve  is  theoretical 
prediction  from  McPhee  (1979). 


remove  inertial  motion.  In.  all,  254  samples  were 
obtained  with  relative  speeds  greater  than  8  cm  s' 1 
(in  order  to  minimize  the  effect  of  uncertainty  in 
F().  Stress  was  calculated  for  each  sample  according 
to  eq  3  with  Cu  s  0.0027  and  m  5  250  g  cm-*.  By 
least-squares  regression,  the  exponent  b  with  a  90% 
confidence  interval,  was  found  to  be  1.78  t  0.12. 
Thus  the  relationship  was  close  to  quadratic  (b  =  2), 
but  statistically  different  at  the  90%  confidence  level. 

The  calculated  magnitudes  of  a  and  the  deflection 
angle  0  between  r0  and  V,  are  sensitive  to  C(  0  and 
m,  both  of  which  have  relatively  large  uncertainties. 
However,  we  showed  (McPhee  1979,  Table  1)  that 
the  exponent  b  is  insensitive  to  variation  in  C, 
and  m  over  a  large  range  of  plausible  values.  We 
therefore  concluded  that,  despite  uncertainties  in 
stress  magnitude  and  0,  the  shape  of  the  r0  vs  V 
curve  as  indicated  by  b  was  well-defined  by  the  data, 
and  that  a  proper  stress  formulation  should  exhibit 
similar  shape.  The  dashed  curve  in  Figure  1  results 
from  the  numerical  model  of  McPhee  (1979):  its 
shape  factor  turned  out  to  be  b  =  1 .70 ±  0.00,  which 
was  within  the  confidence  interval  found  for  the  data. 

It  was  more  difficult  to  assign  an  experimental 
value  to  0.  One  measure  was  the  angle  between  neg¬ 
ative  ice  drift  and  current  velocity  measured  at  0.2 -m 
depth.  Over  the  entire  year,  the  average  for  all 
camps  was  about  24°  (McPhee  1980b,  Table  2)  but 
there  was  much  scatter  in  the  data.  Precise  current 
measurements  from  the  1972  A1DJEX  Pilot  Study 
also  showed  0  to  be  about  24°,  so  this  value  was 
accepted  for  an  average,  but  with  much  less  certainty 
than  the  value  for  b. 

In  the  next  section,  we  consider  several  drag  laws, 
each  with  a  corresponding  simple  PBL  model.  Our 
evaluation  of  each  law  depends  primarily  on  how 
well  it  satisfies  the  shape  requirements  implied  by 
the  AIDJEX  data. 


HIERARCHY  OF  DRAG  LAWS  AND 
SIMPLE  MODELS 

The  drag  law  in  complex  notation 

In  order  to  simplify  as  much  as  possible  a  com¬ 
parison  among  several  drag  formulations,  it  is  useful 
to  express  the  drag  law,  as  follows.  First  define  the 
friction  velocity  as 


where  r  is  the  vector  kinematic  stress  at  the  inter¬ 
faces.  Note  that  ii.  is  a  vector  (complex  number) 
in  the  direction  of  f0  with  magnitude 


u#*u,efa' 


v.ve«a* 


figure  2.  Schematic  of  friction  vehc- 
ity /surface  velocity  relation  in  com¬ 
plex  plane. 


From  comments  in  the  previous  section,  we  ex¬ 
pect  realistic  drag  laws  to  fall  closer  to  quadratic 
(i.e.  tq  «  V2)  than  linear  (rQ  «  V)  and  to  exhibit  a 
turning  angle  of  the  order  of  20°  to  30°.  Using  com¬ 
plex  notation,  we  express  the  drag  law  as  a  ratio  of 
V  and  u »: 


—  =  r  =  r>  <0  . 

u. 


From  Figure  2  we  have 


so  that  r  =  V/u,  is  a  dimensionless  surface  speed  and 
0  =  a2-a,  is  the  rightward  (f  positive)  deflection 
angle  of  surface  velocity  with  respect  to  interfacial 
stress.  Note  that  r  and  0  are  constant  for  quadratic 
drag  with  constant  turning  angle;  all  other  laws  have 
T  and  0  dependent  on  other  parameters  of  the  problem 
(e  g-  «•,/,  etc.).  Expressing  the  drag  law  as  a  com¬ 
plex  ratio  greatly  simplifies  its  algebraic  representa¬ 
tion  and  is  computationally  useful  (e.g.  given  u, , 

V is  found  by  a  scaling,  T,  and  rotation,  e*^). 

The  momentum  equation  for  the 
planetary  boundary  layer 

The  horizontally  homogeneous,  steady-state  PBL 
equation  of  motion  is  given  in  complex  notation  by 

(5) 

where  U{  is  the  velocity  in  a  frame  of  reference  fixed 
to  earth,  t  is  the  kinematic,  turbulent  Reynolds 
stress,  and  we  have  expressed  the  pressure  gradient 
due  to  sea-surface  tilt  in  terms  of  the  geostrophic 
velocity  from  eq  2.  .As  with  the  ice  momentum  equa¬ 
tion,  we  eliminate  V.  by  considering  a  reference 
frame  drifting  with  the  geostrophic  velocity;  hence 
for  V=Ut-  V%,  eq  5  becomes 


3 


if 


(6) 


To  solve  eq  6  an  expression  for  t  is  required  and 
here  we  make  the  simplest  closure  assumption,  i.e. 
that  stress  is  related  to  dll/dz  by  ah  eddy  viscosity 
coefficient: 


bU 
™  02 


Thus  eq  6  becomes 


In  the  remainder  of  this  section  we  show  how  various 
drag  laws  imply  a  variety  of  simple  formulations  for 
the  eddy  viscosity. 


The  last  is  not  necessarily  obvious,  because  turbu¬ 
lence  comprises  a  broad  spectrum  of  length  scales. 
The  interpretation  here  is  that  6  is  the  scale  of  the 
largest  eddies  with  which  smaller  eddies  can  "keep 
up,”  i.e.  diffuse  the  property  to  the  surrounding 
fluid.  Note  that  for  different  properties  in  the  same 
flow,  £  may  have  different  values. 

The  simplest  expression  for  the  flux  of  e  across 
the  horizontal  surface  is  linear  in  each  quantity,  vz 


By  analogy  with  molecular  diffusivity,  we  may  define 
an  “eddy  diffusivity”  Ke  =  v£f  with  units  L2T' 
such  that 


F=-K. 


de 

dz 


Eddy  viscosity 

The  essential  tenets  of  the  eddy  viscosity  concept 
are  that  turbulent  transfer  (overturning  associated 
with  instabilities  in  the  flow  field)  is  much  more 
efficient  than  molecular  transfer,  and  that,  in  some 
way,  this  turbulent  transfer  can  be  expressed  in  terms 
of  mean  properties  of  the  flow.  Note  the  fundamental 
distinction  between  properties  of  the  flow  and  pro¬ 
perties  of  the  fluid,  as  emphasized,  for  example,  by 
Tennekes  and  Lumley  (1972,  Chapter  1). 

Figure  3  depicts  an  overturning  eddy  in  a  fluid 
characterized  by  a  gradient  of  some  property,  e. 

A  conspicuous  feature  of  turbulence  (as  opposed, 
say,  to  organized  waves)  is  its  highly  diffusive  nature; 
fluid  parcels  caught  up  in  a  turbulent  eddy,  as  shown, 
become  involved  in  progressively  smaller  eddies  and 
readily  transfer  their  properties  to  the  surrounding 
fluid.  In  Figure  3  parcel  I  distributes  a  deficit  of  e 
to  the  flow  above  A,  and  parcel  2  distributes  an 
excess  below  A,  so  that  a  downward  flux  of  e  across 
the  surface  z-h  occurs.  This  flux  apparently  de¬ 
pends  on  three  quantities:  first,  the  gradient  of  the 
property  de/dz,  second,  the  overturning  velocity 
v,  and  finally  the  length  scale  £  of  the  most  energetic 
eddies  found  in  that  particular  region  of  the  flow. 


If  the  property  in  question  is  momentum,  then 
the  vertical  momentum  flux,  which  is  minus  the 
horizontal  part  of  the  turbulent  stress  tensor  and 
therefore  a  vector,  is  given  by 

Fm  =  -pT  =  -pAm|f  (8) 

Where  V  is  the  mean  velocity  vector  and  r  =  -  (mV 
+  ivw  )  is  the  Reynolds  stress. 


linear  drag  -  constant  eddy  viscosity 

The  simplest  solution  to  the  PBL  momentum 
equation  is  obtained  when  Km  has  no  depth  depend¬ 
ence  (Ekman  1 905).  With  the  further  stipulation 
that  Km  is  independent  of  other  variables,  the  solu¬ 
tion  provides  a  linear  drag  law.  Eauation  7  becomes 


bHf 

dz2 


with  solution 


U=Ae6z  +&T5* 


i/i 

Wm)  e'"4 


=  <fl2Km)'l*  (1+iV 

As  0  is  bounded  at  depth  (r-*  --),/?  =  0;  therefore 

M(z)=  Ke  e«*  =  Ke  eW2Km),/21 
x  e'W2*nl),/2i 


Figure  3.  Diagram  of  turbulent  flux 
mechanism. 


(9) 


where  VE  is  the  surface  Ekman  velocity. 

Since  5  is  complex,  the  solution  has  both  oscillatory 
and  exponentially  decaying  components  in  z.  pro¬ 
viding  the  well-known  Ekman  spiral  in  which  the 
drift  current  rotates  and  decreays  with  depth.  Note 
that  8  is  the  inverse  of  an  “e-foiding”  PBL  depth 
scale  which  varies  as  ■ 

Our  specific  interest  is  the  relation  between  surface 
stress  and  velocity,  given  by  the  boundary  condition 


-  „  dU 

r0=«.«.=*m 


i=0 


=  (/*„,)' /2  e‘”'A  ■ 

The  drag  law  for  the  Ekman  solution  is  thus 

=  (fKJ-'/iu.e-*”1*  (10) 


-r-  =(A'r)-'/2  u.e-'t5  (11) 

u. 

where  Q  =  20°  and  KR  =  24  cm2  s'1 ,  following  the 
value  for  eddy  viscosity  under  sea  ice  reported  bv 
Hunkins  (1966). 

To  be  consistent  with  Ekman  theory,  reducing 
the  drift  angle  is  tantamount  to  inserting  a  layer 
with  constant  stress  between  the  ice  and  the  Ekman 
layer.  If  this  layer  is  thin  enough,  the  change  in 
stress  across  it  will  be  small,  yet  there  will  be  a  con¬ 
siderable  current  shear  in  the  direction  of  stress.  The 
concept  is  sketched  in  Figure  4,  and  we  shall  analyze 
the  problem  in  some  detail  to  establish  a  pattern  for 
subsequent  refinements. 

V  is  the  vector  sum  of  l'E  and  V%,  the  latter  being 
the  velocity  of  the  surface  with  respect  to  the  average 
velocity  at  the  level  z  =  -h.  Integration  of  eq  8  for 
constant  Km  and  stress  yielus 


i.e.  T  =  U,l(fKm){  f2 , 0  =  45°. 

If  Km  is  constant  with  respect  to  all  variables, 
drawbacks  are  immediately  apparent.  First,  drift 
data  show  that  the  drag  law  is  not  linear,  and  that 
P  is  about  half  45°.  From  the  viewpoint  of  simple 
turbulence  theory,  Km  being  constant  implies  that 
the  product  of  mixing  length  and  turbulent  velocity 
is  a  property  of  the  fluid.  Instead,  we  intuitively 
expect  larger,  more  energetic  eddies  (8  and  v  both 
increasing)  as  the  ice  speed  increases. 

Despite  its  shortcomings,  the  so<alled  Ekman 
approach  has  often  been  used  for  ice  modeling  because 
the  stress  is  linear  in  V.  Some  authors  have  reduced 
the  turning  angle  to  values  closer  to  those  observed. 
For  example,  Rothrock  (1975)  used  a  linear  drag 
law  of  the  form 


where  h  is  the  depth  of  the  constant  stress  layer. 
Note  that  this  complex  ratio  has  no  imaginary  com¬ 
ponent:  Vt  is  in  the  same  direction  as  u..  The 
stress  acting  at  the  top  of  the  Ekman  layer,  z  =  -h, 
is  assumed  unchanged,  so  the  drag  law  is  obtained 
by  combining  eq  10  and  1 1 : 

T  =  —  =(2/Km)-'/2  u. 

u. 

x  [I  -i]  .  (13) 


r«-h 


Constant 


Stress 


Figure  4.  Two-layer  system  with  constant  eddy  viscosity.  Note 
that  shear  in  the  surface  layer  is  constant.  The  actual  current 
rotates  clockwise  with  increasing  depth  in  both  layers. 
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Since  0  is  specified,  the  relation 

tan<J=  1§r)=  (,4) 

determines  h  and  (eq  1 3)  may  be  rewritten  as 

r  =  —  =  (2/Xm)'*'2  csc0u.e-'<>  (15) 

u. 

Note  that  the  eddy  viscosity  »n  the  Ekman  layer  is 
related  to  that  specified  in  eq  11  by 

*m=*Rcsc20/2 

which  for  Rothrock’s  value  (24  cm2  s'1 )  provides  a 
“real”  eddy  viscosity  about  four  times  that  reported 
by  Hunkins.  Solving  eq  14  for  h  with  (3  =  20°  yields 
a  thickness  of  about  10  m  for  the  surface  layer,  which 
is,  of  course,  unrealistically  large  for  the  assumption 
of  constant  stress.  This  exercise,  in  which  the  linear 
drag  law  is  interpreted  in  terms  of  an  Ekman  layer 
overlain  by  a  constant-stress  layer,  demonstrates 
that  such  laws  have  little  basis  in  PBL  theory,  although 
they  are  sometimes  justified  by  Hunkins’  ( 1 966) 
eddy-viscosity  estimate.  It  is  worth  noting  in  this 
regard  that  Ekman  (1905)  treated  eddy  viscosity  as 
constant  only  in  the  sense  of  having  no  vertical  varia¬ 
tion.  In  fact,  he  postulated,  with  characteristic  in¬ 
sight,  that  in  the  open  sea  eddy  viscosity  should  vary 
as  the  square  of  the  surface  wind,  a  view  corroborated 
by  the  present  results. 

Linear  eddy  viscosity  -  the  constant  stress  layer 
Near  a  solid  surface  in  a  developed  turbulent  flow, 
eddy  momentum  flux  is  comparatively  well  under¬ 
stood.  Much  experimental  evidence  indicates  that 
the  proper  velocity  and  length  scales  in  the  eddy 
viscosity  expression  are  v  ~  u,  and  fi  =  Adz  I  where 
k  is  von  Karman's  constant,  k  =  0.4.  Integration 


of  eq  8  for  constant  stress  between  levels  z  =  -h  and 
z  =  -z0  yields  the  familiar  law  of  the  wall: 

V,  =  0(-zo)  -  u(-h) =  \ 111  (”)  "*  ‘  (16) 

The  surface  roughness  length  zQ  represents  a  scale 
at  which  protrusions  on  the  ice  underside  are  sensed 
by  the  turbulent  flow.  For  a  surface  with  non  uniform 
relief  like  an  ice  floe,  it  is  hard  to  assign  a  precise 
physical  meaning  to  z0;  we  take  it  to  be  the  level  at 
which  mean  turbulent  stress  equals  the  overall  inter- 
facial  stress,  and  above  which  there  is  no  mean 
shear  in  the  water  column,  averaged  spatially. 

For  fixed  h,  eq  16  is  a  quadratic  drag  with  no 
turning,  but  the  constant  stress  approximation  and 
the  required  alignment  of  stress  and  relative  current 
clearly  restrict  its  region  of  validity  to  values  of  h 
that  are  small  compared  to  the  total  PBL  depth,  so 
that  V%  is  never  equal  to  V. 

Two-layer  eddy  viscosity 

Figure  5  depicts  a  boundary  layer  in  which  eddy 
viscosity  increases  linearly  until  it  reaches  a  maximum 
at  some  level  z  =  ~h,  below  which  it  is  constant.  If 
the  stress  at  z  =  -h  is  not  much  changed  from  the 
surface  value,  the  water  column  below  h  will  form 
an  Ekman  spiral  as  if  driven  by  the  interfacial  stress. 
Between  the  ice  and  z  =  -h  there  is  a  logarithmic 
velocity  profile  in  the  direction  of  u,.  As  before, 
the  surface  velocity  V  is  the  sum  of  VE  and  K,,  both 
known  functions  of  u.. 

Because  there  is  no  physical  change  in  the  fluid 
at  z  =  -h,  we  assume  that  the  mixing  length  is  con¬ 
tinuous  across  the  transition.  This  still  allows  con¬ 
siderable  freedom  in  prescribing  the  flow.  For  ex¬ 
ample,  we  might  again  specify  that  ATm  is  independent 
of  all  other  variables,  in  which  case  h-  Km/u».  So 
by  combining  eq  10  and  1 6,  we  have 


Figure  5.  Two-layer  system  with  linear  eddy  viscosity  near  surface. 
Shear  in  the  surface  layer  varies  inversely  as  depth,  causing  logarithmic 
current  profile  from  z  *  -h  to  surface. 
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A  problem  is  immediately  apparent  here  because  as 
m.  increases,  the  constant  stress  layer  becomes  thinner. 
The  maximum  size  of  eddies  generated  within  the  sur¬ 
face  layer  (Iz  I  <  A)  thus  decreases  with  increasing 
stress,  so  that  the  drag  law  approaches  the  constant 
Ekman  solution  at  higher  speeds.  All  of  this  is  con¬ 
trary  to  our  physical  reasoning. 

Reed  and  Campbell  (1962)  developed  a  more 
realistic  approach  by  specifying  a  constant  value  for 
A  and  then  matching  the  magnitude  of  the  eddy  vis¬ 
cosity  in  the  Ekman  layer  to  the  surface-layer  value 
at  2  -  -A.  In  effect  they  postulated  a  physical  limi¬ 
tation  on  the  maximum  effective  eddy  size  (mixing 
length)  rather  than  on  the  eddy  viscosity.  Hence 
Km  =  ku,h  where  A  is  an  assigned  constant  and  the 
Reed-Campbell  drag  law  (with  the  simplification 
re  +  A  m  hi)  is 


(2fkh)~112  csc0. 


The  Reed-Campbell  formulation  recognizes  that  eddy 
viscosity  in  the  outer  layer  depends  on  surface  stress, 
but  there  is  still  an  important  conceptual  stumbling 
block:  why  should  A  be  the  same  for  a  flow  with 
low  turbulence  levels  (small  «,)as  for  a  flow  having 
high  turbulence  levels?  We  might  instead  expect  the 
maximum  mixing  length  to  depend  on  other  flow 
parameters  as  well. 

PBL  scaling 

Rossby  (1932)  and  Rossby  and  Montgomery 
(1935)  were  apparently  the  first  to  explore  the  idea 
that  the  maximum  eddy  viscosity  for  the  rotational 
PBL  is  proportional  to  a  velocity  scale  times  a 
length  scale  which  appears  from  dimensional  analysis 
of  the  momentum  equation; namely,  the  velocity 
scale  divided  by  the  Coriolis  parameter /(they  used 
the  geostrophic  wind  speed  in  the  atmosphere, 


analogous  to  V  here).  Much  subsequent  work  for 
the  atmospheric  PBL  has  relied  heavily  on  this  idea 
except  that  the  commonly  used  scales  are  u.  and 
Utlf  for  velocity  and  length  (there  is  a  subtle  differ¬ 
ence  because  u,  and  V  are  not  exactly  proportional). 
The  basic  idea  is  that  turbulent  boundary  layers  in 
rotating  flows  are  dynamically  similar  in  regions 
suitably  far  away  from  solid  surfaces,  provided 
density  is  uniform  and  flow  variables  are  scaled  by 
u.J,  u.,and  u. If  for  kinematic  stress,  velocity, 
and  length  respectively. 

Although  this  framework  has  proved  a  powerful 
tool  for  classifying  and  modeling  the  atmospheric 
PBL,  it  has  not  been  as  widely  accepted  for  oceanic 
modeling  (see,  for  example,  Kraus  1977).  Whatever 
the  reasons  for  this,  the  results  from  under-ice 
studies  of  PBL  mean  flow  and  turbulence  during  the 
AIDJEX  experiments  strongly  support  Rossby-type 
scaling  (McPhee  and  Smith  1976,McPhee  1979  and 
1 980a).  Of  special  significance  are  measurements 
of  turbulent  stress  and  turbulent  energy  spectra  re¬ 
ported  by  McPhee  and  Smith  (1976).  We  found  that 
atmospheric  models  adequately  described  our  tur¬ 
bulence  data  when  scaled  as  above.  Furthermore, 
we  used  peaks  in  the  spectra  of  vertical  velocity 
fluctuations  at  several  levels  through  the  PBL  to  show 
that  a  dimensionless  eddy  viscosity  distribution 
(/Km/u.2)  increased  more  or  less  linearly  with  dis¬ 
tance  from  the  surface  until  it  reached  a  maximum 
at  a  depth  some  fraction  of  the  total  PBL  depth, 
beyond  which  it  decreased  slowly.  While  this  is  not 
the  simple  two-layer  system  discussed  here,  we  found 
(McPhee  1979)  that  there  was  little  difference  between 
numerical  solutions  of  the  PBL  equations  using  an 
exponentially  attenuated  linear  eddy  viscosity  (similar 
to  the  observed  regime)  and  a  linear -constant  two- 
layer  system  (compare  Fig.  1  and  2  of  McPhee  1979). 
This  observation,  in  fact,  prompted  the  present  work. 

Why  turbulent  eddies  should  be  limited  in  size 
by  rotational  effects  (i.e.  that  the  maximum  mixing 
length  should  be  proportional  to  u./f)  is  addressed 
in  some  detail  by  Stern  (1975,  Chapter  8).  From  a 
stability  analysis  of  the  PBL  equation,  he  shows  that 
energy  is  spontaneously  radiated  downward  in  in¬ 
ertial  waves  if  the  scale  of  turbulent  fluctuations 
exceeds  a  length  proportional  to  u,//.  He  concludes 
that  this  provides  an  effective  "brake”  on  surface- 
driven  turbulence  which  tends  to  keep  the  scale  of 
the  largest  disturbances  near  this  instability-governed 
limit. 

A  dimension  less  two-layer  system 

The  ideas  above  can  be  easily  incorporated  into 
the  two-layer  model  as  sketched  in  Figure  6.  All 


Figure  6.  Two-layer  system  nondimensionalteed  by  Rossby-similarity 
scales.  Note  that  actual  eddy  viscosity  in  Ekman  layer  is  proportional 
to  u.2/f. 


quantities  are  nondimensionalized  and  we  assume 
that  there  is  a  “universal”  maximum  nondimensionaJ 
mixing  length  {N.  The  outer  layer  eddy  viscosity 
is  then 


and  the  Ekman  solution  (eq  10)  is 

—  -  (**N)',/J  e  '*/4 

it. 

which  is  clearly  a  quadratic  drag  law.  The  surface 
layer  solution  is 


K 

it. 

where  *fzju.  is  the  “nondimensional  roughness 
length.”  The  drag  law  for  the  combined  system  is 


r-m»)r{l'ctct.  (21) 

Equation  19  reduce*  to  a  pure  quadratic  drag  law 
if  t0  it  constant.  This  stipulation  is  questionable, 


1  ./In' 


irln  ^ 


because  unlike  the  turbulence  scales  in  the  outer 
layer,  z0  is  a  physical  property  of  the  surface  itself, 
and  there  is  no  a  priori  reason  to  assume  it  scales 
with  u.lf.  In  formulating  the  drag  law  for  the 
AIDJEX  Model  (McPhee  1975),  we  assumed  that 
over  the  limited  range  of  u.lf  encountered  during 
meaningful  simulations,  the  effect  of  variations  in 
the  “surface-friction  Rossby  number”  (Ro.  =  u.lfza 
=  1  /£0)  would  be  minor  and  that  a  quadratic  drag 
law  would  suffice.  This  was  in  part  due  to  uncertainty 
about  a  suitable  value  for  za,  since  it  represents  an 
aggregate  roughness  over  perhaps  several  large  floes. 
The  success  of  subsequent  simulation  of  summer 
drift  using  the  quadratic  drag  (McRiee  1980b)  tended 
to  reinforce  this  view. 

On  the  other  hand,  analysis  of  drift  statistics 
described  previously  indicated  that  the  drag  law  is 
slightly  different  from  quadratic,  so  we  might  rather 
treat  z„  as  constant,  which  gives  both  0  and  V  slight 
dependence  on  u.lf  through  the  term  In  (£N/£0) 
in  eq  20.  The  two-layer  system  of  eq  19  with  z„ 
constant  includes  most  of  the  features  required  for 
a  conceptual  model  in  which  die  size  of  momentum 
exdunging  eddies  increases  as  the  distance  from  the 
interface,  but  is  limited  to  some  fraction  of  u.lf  by 
an  instability  mechanism  which  rapidly  propagates 
energy  away  from  the  boundary  layer  once  the  size 
limit  is  exceeded.  All  such  flows,  when  scaled  as 
above,  are  similar  except  in  a  thin  region  right  next 
to  the  interface  in  which  the  other  physical  length 
scale  z0  asserts  itself. 

As  an  historical  note,  Ekman  (1905)  assumed  that 
the  depth  of  frictional  influence  [D  ■  *(2 Kmlf)112) 
in  the  open  ocean  varied  as  W/( sin  9)'/*  where  W 
is  the  surface  wind  velocity  and  9  is  latitude.  We 
see  that,  if  the  dependence  were  rather  D  «  W/ sin  9 
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(proportional  to  ujf),  this  would  provide  an  eddy 
viscosity  such  that  fKmlu»2  is  constant  (rather  than 
JCm),  resulting  in  a  quadratic  drag.  If  the  minor 
variation  of /at  high  latitudes  is  ignored,  faithful 
application  of  Ekman’s  ideas  does  result  in  a  quad¬ 
ratic  drag,  and  the  turning  angle  becomes  realistic 
when  a  surface  layer  such  as  that  of  Shuleikin  (1938) 
is  provided. 

While  the  modeling  results  of  McPhee  (1979) 
showed  that  the  particular  shape  of  the  nondimensional 
eddy  viscosity  distribution  in  the  outer  layer  is  not 
overly  important,  there  remains  an  obvious  question 
about  the  two-layer  approach:  specifically ,  how 
valid  is  the  assumption  of  constant  stress  in  the  sur¬ 
face  layer?  Here  the  numerical  results  are  less  en¬ 
couraging:  they  show  that  it  the  level  £  =  -fN,  the 
nondimensional  stress  has  decreased  in  magnitude 
by  about  20%  and  rotated  clockwise  perhaps  20° 
from  its  surface  value  (sec  Fig.  2  of  McPhee  1 979). 

This  suggests  a  final  modification  to  the  two-layer 
system. 

A  dimensionless  two-layer  system  with  modified  stress 
Equation  6,  nondimensionalized,  has  the  form 

*  =  |  (22) 

where  £  =/z/«.,u  =  U/u and  T*  r/u.u..  The  non- 
dimensional  stress  at  $  is  thus 

o 

7tt)=ro-fJ  udi  (23) 

t 

where  the  subscript  o  denotes  the  surface  values.  To 
find  jE  ■  7(-fN  ),  we  first  note  that  T  is  aligned  at 
the  surface  with  the  real  axis  (i.e.  f„  *  1),  and  then 


assume  that  within  the  surface  layer  (1$  I  <  $N)  shear 
in  the  lateral  (imaginary)  component  of  velocity  is 
negligible.  Since  the  real  part  of  the  surface-layer 
velocity  varies  logarithmically,  integration  of  eq  23, 
evaluated  with  the  approximation  £N  +  £„  *  £N , 
yields 

rE*l -i(vE  +  l/*)fN  (24) 

where  vE  =  m(-|n)  is  the  nondimensional  Ekman 
velocity  forced  by  the  attenuated  and  rotated  surface 
stress  TE.  Using  this  as  the  upper  boundary  conu- 
tion  in  the  nondimensional  Ekman  solution  we  have 

vE  =  ir  =  (*«n)1/2  (25) 

«• 

and  substituting  into  eq  24 

rE=(i-ftN/*)/li+(«N/2*),/J(i  +  01 

*  1  -(iN/2*)1/2  (1  +0  -  (26) 

For  the  value  £N  =  0.045,  derived  in  the  next  sub¬ 
section,  we  have  7*E  =  0.80, 0  =  -17.3°,  which  agrees 
well  with  the  numerical  result  of  McPhee  (1979). 
Solving  for  vE  from  eq  24,  we  have 

VE  =(2*6n)»/2  (I -o-  1/* 

r~  ~  ~ (2*$n)_,/2  {i+(2 fN/*)WI 

»• 

X  NWU-'I  -*■}  <27> 

and  the  system  is  sketched  in  Figure  7. 


C  hC, 


Figure  7.  Two-kyer  system  in  which  street  is  allowed  to  my 
through  a  logarithmic  surface  kyer.  The  direction  end  rektive 
magnitude  of  the  friction  velocity  at  the  top  of  the  Ekman  kyer 
it  Moated  by  ii.E.  Numbers  betide  current  vectors  denote 
nondimensional  depth  units,  -fa/u,. 


Figure  8.  Nondimensional  stress  and  velocity  in  the  planetary 
boundary  layer  according  to  two-layer  model  with  modified 
stress  (solid  curves)  and  atmospheric  PBL  model  of  Businger 
and  Arya  (1974)  (dashed  curves).  Velocity  provile  in  the  sur¬ 
face  layer  (\$N\  <  0.045)  depends  on  Zq. 


In  addition  to  the  drag  law,  this  approach  allows 
a  complete  specification  of  velocity  and  stress 
throughout  the  boundary  layer.  The  E  Ion  an  solution 
for  the  outer  layer  is 

M<«)=vEe’,(tMN),  *<-«N  (28) 

where  rj  =  u.8//=  (I  + i).  For  the  surface 

layer, 

=  ln(iN/l«l>. 

lo«lll<*M-  (29> 

The  nondimensional  stress,  T  =  r/u.u.,  is,  for  the 
outer  layer, 

7tt)=rEe,(,  +  *N\  *<-{N  (30) 

and  for  the  surface  layer  the  stress  is  well  approxi¬ 
mated  by 

itt)-l+l(»E  ♦!/**. 

«o«l*l<*N-  CO 

Equations  28-31  constitute  a  closed  PBL  model  for 
the  neutral  boundary  layer  unaffected  by  a  pycnocline. 
Solutions  for  {  <  -(N  (=  0.04S)  are  drawn  in  Figure 
8  along  with  results  for  the  neutrally  stable  atmos¬ 
pheric  PBL  from  the  numerical  model  of  Businger 
and  Arya  (1974).  The  agreement  is  remarkable  con¬ 


sidering  that  the  only  adjustable  parameter  in  the 
analytic  solution,  (;N,  was  chosen  by  comparing  the 
drift  of  sea  ice  with  the  stress  exerted  upon  it  by  the 
wind. 

Evaluating  the  drag  laws 

From  the  concept  of  a  surface  layer  separating 
the  ice-ocean  interface  from  a  pure  Ekman  layer, 
we  have  constructed  a  sequence  of  six  drag  laws, 
each  one  hopefully  including  a  slightly  more  sophis¬ 
ticated  view  of  the  boundary-layer  structure  than 
the  one  before  it.  Most  expressions  for  water  stress 
used  in  ice  modeling  fit  into  one  of  these  categories: 

1.  Constant  eddy  viscosity,  reduced  turning 
angle  (eq  13). 

2.  Logarithmic,  constant-stress  surface  layer, 
constant  eddy  viscosity  (eq  1 7). 

3.  Reed-Campbeli  approach,  h  constant  (eq  18). 

4.  Nondimensional  eddy  viscosity,  constant  fi 
(eq  19,  constant). 

5.  Rossby  similarity  scaling,  constant-stress 
surface  layer  (eq  19,  z0  constant. 

6.  Rossby  similarity  scaling,  logarithmic  surface 
layer  with  modified  stress  (eq  27). 

In  order  to  compare  the  drag  laws  among  them¬ 
selves  and  with  the  AIDJEX  free-drift  data,  values 
for  the  pertinent  parameters  in  each  expression 
were  chosen  so  that  when  u*  is  equal  to  1  cm  $'* , 

T  *  1 3  and  fi  *  24°.  The  value  for  T  comes  from 
the  fitted  r  vs  V  curve  (Fig.  1)  evaluated  at  u*  =  1 
cm  i-1 ,  which  was  close  to  the  average  condition 
encountered  during  the  1975  melt  season.  Ex¬ 
pressions  for  fi  and  T  are  summarized  in  Table  1 


Table  I .  Summary  of  drag  laws. 


Drat  law* 

r 

Parameters 

Stress-speed  equation 

1 

24° 

(2/Xm)-*/a  csc/3u* 

Km  =  128  cm1  s'1 

t  =  0.077  V 

2  cot'1 

•  +W.),,a  ,  / 

I+  ku. 

fe)  (W^csc*. 

X_  =  128  cm1  s'1 
m 

r  =  0.023  Ff‘-46±00J) 

z0  =  23  cm 

3  cot'1  | 

,  JW\112  i_  h 

Zo 

|  (2/JWi)",/J  c&cfju*1!2 

h  =  319  cm 

T  =  0.0177F<,s7tOOO> 

z0  =  23  cm 

4 

24® 

13 

r  =  0.0059  F7 

5  cot'1  11  + 


6  cor1 


1  + 


W"-r)| 


(2tfN)-‘'J  csc/J  {N=  0.045  r  =  0.0128l*,-7°to-oo> 

zQ  =  23  cm 

(2^N)-«/2  csc/3  fN  =0.045  r  =  O.0128^,  7«»±0  00) 


zQ  =  8.5  m 


~  *  Te-'4 
U* 


for  each  model.  Note  that  no  attempt  was  made 
to  choose  “physically  realistic*’  parameters  garnered 
from  the  literature.  Our  reasoning  was  that  the 
laws  were  best  compared  by  having  them  all  agree 
at  one  point  within  the  range  of  acceptable  data. 

Curves  coinciding  with  the  chosen  parameters 
are  plotted  in  Figure  9  along  with  the  envelope  of 
the  experimentally  determined  exponent,  b  *  1 .78 
i  0.12.  In  addition,  these  curves  were  fitted  to  the 
stress  relation 

rmaVb 

with  a  least-squares  analysis  of  the  29  points  used 
to  generate  each  curve .  Values  of  b  with  90%  con¬ 
fidence  brackets  are  tabulated  for  each  of  the  drag 
laws  in  the  last  column  of  Table  1 .  Only  laws  5 
and  6  fall  within  the  90%  error  bounds  of  the  ob¬ 
served  data.  The  stress-speed  relationship  is  thus 


best  described  by  the  Rossby-similarity  laws,  which 
include  the  most  extensive  treatment  of  the  prin¬ 
ciples  governing  turbulent  flow.  Note  that,  although 
the  drag  results  are  very  similar  for  both  Rossby- 
similarity  approaches,  the  effective  surface  rough¬ 
ness  is  smaller  in  the  reduced  stress  model  (eq  27), 
falling  close  to  the  value  (10  cm)  determined  in 
McPhee  (1979).  As  we  have  shown,  the  two-layer 
model  with  reduced  stress  is  a  close  approximation 
of  the  numerical  results  of  Businger  and  Ay ra  (1974) 
and  McPhee  (1979)  for  determining  stress  and  velocity 
at  all  levels  in  the  boundary  layer.  In  flows  where 
a  shallow  pycnodine  occurs,  flow  profiles  are 
affected  at  lower  levels  and  drag  characteristics 
are  apparently  changed  slightly  as  discussed  in  McPhee 
(1979).  We  reiterate  that  when  the  ice  pack  is  weak 
enough  to  admit  inertial  oscillation,  the  drag  laws 
here  are  appropriate  only  for  motion  averaged  for 
time  scales  of  the  same  order  or  longer  than  the 
inertial  period  (~  1 2  hours). 
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Figure  9.  Comparison  of  six  drag  law  formations  with  data  from 
1975  AIDJEX  Experiment.  Shaded  portion  shows  the  limits  of 
the  empirical  fit  included  in  the  90%  confidence  interval. 


Rossby  similarity  parameters  and  buoyancy  effects 

The  usual  approach  to  parameterizing  surface 
drag  on  the  atmospheric  free  stream  using  Rossby- 
similarity  scaling  is  somewhat  different  from  that 
developed  in  the  previous  section.  Blackadar  and 
Tennekes  (1968)  summarize  the  technique,  which 
matches  the  mean  shear  in  an  overlap  region  where 
the  widely  different  scales  z„  and  u.//both  affect 
the  flow.  Since  the  nondimensional  shear  is  the 
same  by  either  scaling,  a  logarithmic  velocity  pro¬ 
file  ensues,  which  is  matched  to  the  empirical  surface- 
layer  profile.  The  result  is  usually  expressed  in  terms 
of  the  nondimensional  velocity  components  parallel 
and  perpendicular  to  the  surface  stress,  i.e. 

Re  (r)  =  Re  ^  =  l-  (In  /to.  -  A )  (32) 

lm(0=-f-  (33) 

where  /to.  =  1/|0  =  u.//z0,  as  before,  and  A  and  B 
are  found  empirically.  The  last  two  drag  laws  (eq 
19  and  27)  give  the  same  result  as  the  Blackardar- 
Tennekes  approach:  e.g.  for  eq  27 

Ifc  (h-(2*tN )-'/»+ £(ln  In -o 

+  j-  InJto. 

imav^iN)-1'*. 


Hence 

*  =  (*/2|n),/j 

A  =  -(B  +  ln  |N-1). 

For  In  =  0.045,  the  parameters  are  (for  neutral 
stability)  A  *  1 .99, B  -  2.1 1 .  Conversely,  the  drag 
parameters  T  and  0  may  be  expressed  in  terms  of 
A  and  B  by 

0=tan->  j*/(ln/to.  -A)\  (34) 

r  =  £  { B*  +(ln/to. -A)2)'!*  .  (35) 

In  naturally  occurring  boundary  layers,  buoyancy 
forces  associated  with  turbulent  fluctuations  in  a 
nonuniform  density  field  often  play  an  important 
role  in  flow  dynamics.  Generally  speaking,  in  a 
horizontally  homogeneous  flow,  two  new  length 
scales  (in  addition  to  z0  and  «.//)  become  important 
when  buoyancy  effects  are  considered.  They  are 
the  depth  of  the  mixed  layer  and  the  Obukhov  length 
(see  Arya  and  Sundararajan  1976,  for  a  discussion 
of  the  atmospheric  analogs).  While  the  former  is 
certainly  important  when  it  is  small  compared  to 
u.//,  it  does  not  fit  easily  into  the  simple  theory 
here.  Our  previous  modeling  results  (McPhee  1979) 
suggest  that  the  effect  of  the  pycnocline  on  mean 
surface  drag  is  not  large  under  commonly  encountered 
conditions. 


y,*5 
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The  second  length  scale,  L  (the  Obukhov  length), 
is  used  to  quantify  the  effect  of  surface  buoyancy 
flux  due  to  ice  growth  or  ablation.  L  is  obtained 
(see,  for  example,  Tennekes  and  Lumley  1972, 
Chapter  3)  by  nondimensionalizing  the  surface-layer 
turbulent  kinetic  energy  equation  such  that  the  non- 
dimensional  production  rate  of  turbulent  kinetic 
energy  by  buoyancy  at  a  particular  level  z  is  given 
by  z/L  where 

L  =,P0u,3l(gk  pV)  (36) 

with  pV  being  the  turbulent  flux  of  density  fluc¬ 
tuations.  When  L  is  positive,  turbulence  levels  are 
reduced  because  the  eddies  must  work  against 
gravity ;  when  L  is  negative  turbulence  is  enhanced. 
The  magnitude  of  L  indicates  the  scale  at  which 
buoyancy  is  important:  e.g.,  if  |  L  |  is  very  large 
(pw  small)  buoyancy  has  little  effect  on  flow 
scaling.  On  the  other  hand,  if  |£|  is  much  less  than 
the  scale  that  turbulence  would  have  in  the  absence 
of  density  flux,  then  the  turbulent  structure  may 
be  completely  dominated  by  buoyancy. 

Without  examining  the  complex  details  of  buoy¬ 
ancy-dominated  flows,  we  shall  make  some  general¬ 
izations  relating  to  freezing  rates  based  on  scaling 
and  show  how  a  simple  extension  of  the  theory  in 
the  last  section  has  application  to  conditions  of 
rapid  melting. 

The  equation  of  state  for  sea  water  is  approxi¬ 
mately 


When  ice  melts  (d< 0)  due  to  contact  with  warm 
water,  the  interfacial  heat  flux  is  pQcp Ad  where  A 
is  the  kinematic  latent  heat  of  fusion.  Thus 

'wW  I  =  -Ad  (</<0|.  (37) 

O 

The  nondimensional  Obukhov  length,/..  =/l/u.,  is 

/,.=£- HI A,  +  aA)'*^j-.  (38) 

For  temperatures  less  than  0°C,  the  thermal  expan¬ 
sion  coefficient  a  is  negligibly  small  (but  always 
positive  for  S  >  24°/oo);  for  water  at  6  -  1°C, 

S  *  30°/oo,  or  a  5x  10_s  K~l  (Neumann  and  Pierson 
1966). 

Using  typical  values  (fi-  0.8,  pice/p0  =  0.9, 
s»  -^ce  =  25°/oo,  A  =  65  K,  /=  1 .4  x  10~«  s'1) 
in  eq  38,  for  0W  <  0°C 

L,  =  (-2.0  x  10-s  cm-1  s)  •  u}!d  .  (39) 

For  0W  -  1°C,  the  magnitude  of  L,  is  increased  by 
about  20%  by  the  thermal  expansion  term. 

From  the  scaling  considerations  discussed  above, 
we  expect  L.  to  have  an  impact  on  the  turbulent 
regime  when  its  magnitude  is  of  the  same  order  or 
smaller  than  the  neutral  nondimensional  PBL  depth, 
often  taken  to  be  about  O.S.  For  growth  rates  such 
that 

|cf|  >  (4  x  10"s  cm-1  s)  u.2 


*-«(®-0o )  +  *S-S„) 

Po 

where  p0  is  a  reference  density  at  temperature  0o 
and  salinity  S0;  hence. 


*  -a  wV  +  0  wV 
P  o 

where  wV  and  wV  are  turbulent  fluxes  of  temper¬ 
ature  and  salinity.  When  sea  water  freezes,  only  a 
small  portion  of  the  dissolved  solids  (typically  about 
1/6)  are  retained  in  the  sea  ice  matrix,  so  that  a 
salinity  flux  occurs  as  given  by 


wy  >#  p0  1000 


where  rf  is  the  ice  growth  rate  and  5W  -  Jlc,  is  the 
salinity  difference  in  per  milk.  When  ice  is  rowing 
(d>0),  heat  exchange  is  entirely  latent  and  w’f  m0. 


buoyancy,  therefore,  may  affect  surface  drag.  For 
typical  conditions  («.  ~  1  cm  s-1)  this  effect  is 
about  4  cm  per  day.  But  freezing  rates  this  high 
occur  only  at  low  temperatures  in  comparatively 
thin  ice.  Under  normal  circumstances,  the  drag 
laws  discussed  here  apply  to  pack  ice  that  is  fairly 
thick  (at  least  1  m)  and  compact  (say  80%  cover), 
so  an  order-of-magnitude  estimate  suggests  that  in¬ 
stability  associated  with  surface  freezing  will  rarely 
have  much  effect  on  the  average  ice  drag.  Similarly, 
melting  rates  encountered  in  the  perennial  pack  of 
the  central  Arctic,  where  oceanic  heat  flux  is  small, 
cannot  be  consistently  as  large  as  4  cm  day'1 ,  or 
most  of  the  ice  would  disappear  over  a  60-day  melt 
season. 

On  the  other  hand,  there  are  conditions  under 
which  very  rapid  melting  might  have  a  decided  sta¬ 
bilizing  effect  on  the  PBL  An  example  is  advection 
by  wind  of  pack  ice  across  frontal  zones  at  the  ice 
mar^n  where  mixed-kyer  temperature  may  change 
by  as  much  as  2  to  3  K  within  several  kilometers. 
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Figure  10.  The  Rossby-simikrity  parameters  A  and  B  as 
functions  of  p,  =  1/L„  Data  points  with  error  bars  are  from 
Clarke  and  Hess  (1974 );  open  circles  are  from  model  of 
Businger  and  Arya  (1 9  74);  curves  are  from  present  theory 
with  (Rf)max  «  0.2  (solid)  and  (Rf)max  =  0.15  (dashed). 


There  is  some  indication  at  the  Bering  Sea  ice  margin 
that  melting  rates  at  least  an  order  of  magnitude 
greater  than  4  cm  day*1  occur  when  ice  moves 
across  such  fronts*.  If,  for  example,  we  were  inter¬ 
ested  in  estimating  the  distance  that  thick  floes 
might  migrate  into  warm  water  before  melting,  then 
buoyancy  would  have  to  be  considered  in  calculating 
drag.  The  remainder  of  this  section  addresses  this 
problem  as  an  extension  of  the  mixing  length  approach 
developed  above. 

The  nondimensional,  two-layer  theory  for  a 
stably  stratified  flow  follows  from  the  idea  that  the 
maximum  nondimensional  mixing  length  is  fjN  when 
L,  is  infinite  (neutral  stability),  but  is  proportional 
to  L,  under  very  stable  conditions  (L.  small).  The 
simplest  expression  possessing  these  limits  is 

=  0  *  (40) 

where  X  is  a  constant  yet  to  be  determined.  Given 
L.,  we  sutstitute  £m  for  fN  in  eq  26  and  27  to 
determine  the  stress  attenuation  and  to  solve  for 
the  drag  parameters  T  and  0,  from  which  the  velocity 
may  be  obtained. 

To  evaluate  X,  we  follow  the  reasoning  of  Zilitin- 
kevich  (1975).  With  the  surface-layer  approximation 


*Ptnonal  communication  with  C.  Peat*.  NOAA  Pacific 
Marina  Environmental  Laboratory,  1981. 


(f  =  u.u«),  the  flux  Richardson  number  (i.e.  the 
ratio  of  buoyant  production  to  stress  production  of 
turbulent  kinetic  energy  -  see  Tennekes  and  Lumley 
1972,  p.  98)  is  given  by 

«<--(£  trKfe)  h- 

The  nondimensional  eddy  viscosity,  ( fK(um  J)  in  the 
surface  layer  is  so  we  have  £  *  RfL,.  Observations 
show  that  Rf  has  an  upper  limit  of  about  0.2.  For 
L,  small,  eq  40  approaches 

(•m  at  XL. 

indicating  that  X  =  (^f)„,,x  *  0.2. 

A  comparison  of  the  present  theory  for  two  dif¬ 
ferent  values  of  X  with  other  results  is  shown  in 
Figure  10,  where  the  Rossby  similarity  parameters, 
as  defined  by  eq  32  and  33,  are  plotted  as  func¬ 
tions  of  the  variable  p.  =  1  \L%.  Also  shown  are 
predictions  of  Businger  and  Arya  (1974),  whose 
model  uses  a  nondimensional  eddy  viscosity  based 
on  the  log-linear  wind  profile  determined  from  the 
Kansas  surface-layer  experiment  (Businger  et  al. 

1971). 

The  data  points  with  error  bars  in  Figure  10  are 
from  an  analysis  by  Clarke  and  Hess  (1974)  of  wind 
profiles  throughout  the  entire  PBL  under  stable 


Figure  11.  Data  from  Clarke  and  Hess  (1974)  after 
adjusting  u ,  by  a  constant  factor  to  make  neutral 
values  for  B  agree,  compared  with  present  theory  with 
-  0-2- 


conditions  caused  by  radiative  cooling  during  the 
Wangara  experiment  in  the  Australian  desert. 

By  applying  the  same  arguments  used  to  deduce 
the  neutral  boundary -layer  model  to  a  regime  in 
which  the  maximum  mixing  length  is  given  by  eq 
40,  it  can  be  shown  that  the  curves  in  Figure  10  are 
given  exactly  by 

*-(*/2*,,),,20+W.),/a  (41) 

A  =  l-ln|N-(A:/2|N)1/1  (1  +  7M.)1/2 

+  In  (1  +  7/r.)  (42) 

where  p.  =  l/Z..  and  7  =  $N/(*fWx  • 

For  £n  =  0.045  and  (/?f)mtx  =  0.2,  they  are 

£=2.11 (1  +0.22p.),/2  (43) 

A  =4.10 -2.11  (1  +  0.22  p.)' 

+  In  (1  +0.22  p*)  .  (44) 

Although  at  first  glance,  results  for  B  shown  in 
Figure  10  might  indicate  that  (£f)mlx  =  0.2  is  too 
large,  such  an  interpretation  must  be  approached 
with  caution.  Clarke  and  Hess  (1974)  point  out 
that  their  data  derive  from  scaling  by  u.  as  de¬ 
termined  at  a  very  smooth  site,  and  one  may  infer 
that  the  u,  to  which  the  entire  PBL  responds  is 
larger.  This  point  is  illustrated  in  Figure  1 1 ,  where 
we  have  adjusted  their  data  values  by  applying  a 
constant  correction  factor  to  u*  such  that  their 
neutral  value  (B  *  4.3)  agrees  with  ours  (B  *  2.1 ). 
We  have  assumed  that  geostrophic  wind  and  surface 
heat  flux  measurements  are  representative,  so  that 
ordinate  and  abscissa  values  are  divided  by  about  2 


and  4,  respectively.  Apparently,  the  data  are  too 
scattered  for  one  to  determine  appropriate  values 
for  (£f)mlx.  but  they  do  exhibit  the  same  general 
form  as  eq  41  and  42. 

To  recap,  given  interfacial  stress  and  ice  growth 
rate,  we  can  determine  L.  from  eq  38  and  use  the 
similarity  expressions  eq  43  and  44  to  find  /?  and  T 
from  eq  34  and  35.  For  a  concrete  example,  let  u * 

=  1  cm  s"1 ,  d  *  50  cm  day* ,  zD  =  9  cm,  and /=  1 .23 
x  10~4  s-1  (latitude  58°N);  then  L.  =  0.037,  B  =  5.6, 

A  -  0.45,  T  =  21. 2,  and  0  =  41°.  Compared  with  the 
neutral  case,  the  ice  travels  more  than  half  again  as 
far  under  similar  stress  conditions  because  of  the 
rapid  melting. 

There  is  a  final  ad  hoc  addition  to  the  theory  for 
which  there  is  at  present  little  experimental  justifica¬ 
tion,  but  which  may  be  helpful  in  interpreting  data 
on  drift  and  melt  rates  near  the  margins  of  ice  packs. 

It  is  useful  to  consider  the  physical  scales  implied  in 
the  specific  problem  above.  Regardless  of  how  the  PBL 
depth  is  defined,  similarity  of  profiles  in  the  outer 
(Ekman)  portion  requires  that  it  be  proportional  to 
(kfm)'  I2 .  (Incidentally,  this  implies  that  the  PBL 
depth  is  proportional  to  p»~l  1 2  for  small  Z..,  which 
is  again  consistent  with  Businger  and  Arya  (1974) 
and  Zilitinkevich  [1975]).  In  the  example  above, 
if  the  neutral  PBL  depth  is  0.5  u.//=  41  m,  then 
the  stratified  PBL  depth  is 


=  15.4  m 


so  that  the  thickness  of  the  boundary  layer  is  much 
reduced.  We  may  ignore  the  surface  logarithmic 
layer,  since  $m  u, If  is  only  about  40  cm  (it  varies 
asji*-1  for  large  p,). 
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There  is  no  steady-state  thermal  boundary  layer 
formulation  that  corresponds  with  that  of  the  tur¬ 
bulent  momentum  flux  in  the  outer  layer;  neverthe¬ 
less,  the  maximum  mixing  length  for  temperature 
fluctuations  should  be  of  the  same  order  as  for 
momentum,  and  we  have 


w'd' 


u'2k*m  b± 

f  dz  ■ 


If  we  assume  further  that  the  temperature  gradient 
is  proportional  to  A0IDi  where  AO  is  the  difference 
between  the  ice  melting  temperature  and  the  temper¬ 
ature  at  depth  Dt,  we  obtain  an  explicit  relationship 
between  AO  and  melting  rate,  i.e. 


</  =  - 


k$m  (u.2lf)  AO 
C|  AD$ 


^  C|*$N 
A 


where  we  have  used 


(0.5  u.lf). 


The  constant  c,  is  unknown,  but  ought  to  be  of 
order  unity,  in  which  case 


For  the  example  above,  AO  =  -2.8  K,  which  is  not 
unreasonable  for  conditions  sometimes  found  near 
the  ice  edge.  Note  that  the  depth  of  a  pre-existing 
mixed  layer  is  irrelevant  unless  it  happens  to  be 
less  than  Dt ,  assuming  that  the  ice  advects  directly 
into  a  mixed  layer  of  uniform  temperature  above 
the  melting  point. 

Within  the  theoretical  framework,  the  constant 
c,  could  be  evaluated  experimentally  by  observing 
sea  surface  temperature  AO,  surface  wind  stress 
(u.)and  ablation  rates  of  large  floes  drifting  across 
the  ice  edge  front.  Whether  conditions  steady  or 
horizontally  homogeneous  enough  to  provide  suit¬ 
able  measurements  can  actually  be  found  is  un¬ 
certain;  but  if  C|  were  known,  knowledge  of  wind 
conditions  and  temperature  change  across  the  front- 
both  of  which  might  be  sensed  remotely-would  be 
enough  to  make  a  quantitative  estimate  of  how  long 
floes  would  survive,  how  far  they  would  drift,  and 
to  what  extent  they  would  cool  the  ocean. 


DBCUSSfON 


Rossby-simiiarity  scaling  are  consistent  with  the 
averaged  stress  magnitude  vs  relative  ice  speed  rela¬ 
tionship  determined  from  the  summer  AIDJEX  data. 
When  such  scaling  is  applied  to  the  classical  two-layer 
approach  to  the  PBL  problem,  with  the  added  feature 
that  stress  is  allowed  to  vary  realistically  through  the 
inner  (surface)  layer,  an  analytic  solution  is  found 
that  reproduces  closely  the  numerical  results  of 
McHiee  (1979)  and  of  Businger  and  Arya  (1974). 

The  eddy  viscosity  in  the  outer  (Ekman)  layer  is 
given  byAm=A:{Nu.2//with  {N,the  maximum 
dimensionless  mixing  length,  being  about  0.045. 

With  this  value,  the  drag  law  we  propose  is 

r  =  (p  =2.5  In  Ro.  -4.98)  -5.27/  (45) 

where  Ro,  =fzju,.  The  effective  zQ  inferred  from 
the  AIDJEX  results  is  about  9  cm,  but  may  vary  de¬ 
pending  on  ice  characteristics  such  as  ridging  intensity 
or  previous  mel'ing. 

When  drifting  ice  is  also  melting  rapidly,  buoyancy 
affects  the  drag.  This  tendency  is  modeled  by  a 
simple  extension  of  the  two-layer  theory  in  which 
the  maximum  mixing  length  is  reduced  from  its 
neutral  value  as  stability  increases,  eventually  becoming 
proportional  to  the  Obukhov  length.  With  the  pro¬ 
portionality  constant  equal  to  (/?f)mIX  =  0.2,  the 
results  are  similar  to  those  of  Businger  and  Arya 
(1974)  and  are  supported  by  measurements  of  geo- 
strophic  similarity  parameters  (A  and  B)  in  the 
stable  atmospheric  PBL  reported  by  Clarke  and  Hess 
(1974).  Our  suggested  treatment  of  the  stable  PBL, 
given  the  growth  rate  d,  is 

L.  =(-0.17cm-1  s  2)(/u.*ld) 

5  =  2.11  (1  +0.22 ./!.)»/* 

A  =4.10-2.11  (1  +  0.22/L,)112 
+  In  (1  +0.22  /L.) 

=  jr 4(ln /?o. -A)-iB\  . 

u.  * 

The  results  of  a  specific  example,  3,  which  showed 
that  rapid  melting  might  increase  the  ice  speed  by 
half  its  neutral  value  under  similar  stress  conditions, 
allows  us  to  speculate  on  the  origin  of  an  interesting 
phenomenon  observed  near  the  ice  margin  in  the 
Bering  Sea.  According  to  C.  Pease  and  R.  Bourke* 
it  is  not  uncommon  to  encounter  a  continuous  band 
of  jumbled,  compact  floes  parallel  to  the  edge  of  the 


An  evaluation  of  several  sea -ice  drag  formulations 
shows  that  only  those  conforming  in  some  way  with 


'Personal  communication  with  C.  Paaaa,  NOAA  PM  EL,  I9SI 
and  R.  Boarka,  Naval  Poctaradnata  School.  |9S0. 


main  pack,  but  separated  from  it  by  several  kilo¬ 
meters  of  open  water.  Satellite  imagery  also  shows 
such  features  occasionally  (e.g.,  see  Fig.  3  of  Muench 
and  Charnell  1977).  Results  given  in  the  previous 
section  suggest  the  following  scenario  for  the  for¬ 
mation  of  such  features.  Prior  to  the  onset  of  an 
off-ice  wind  event,  the  ice  at  the  margin  of  the  pack 
is  consolidated  and  ridged,  apparently  by  wave  action 
and  other  edge  effects.  As  the  wind  begins  blowing 
toward  Open  water,  the  entire  pack  moves,  but  the 
leading  floes,  drifting  into  warmer  water,  accelerate 
away  from  the  main  pack,  cooling  the  water  behind 
them.  Thus  the  leading  edge,  which  was  thicker  to 
begin  with,  experiences  less  drag  at  all  times,  and 
separates  from  the  main  pack  more  or  less  intact. 

As  the  wind  dies,  heat  flux  (which  depends  on  u.) 
and  melting  decrease,  so  that  the  band  may  then 
persist  for  some  time  as  a  stationary  feature. 

Whether  or  not  such  a  situation  occurs,  the 
framework  presented  here,  which  has  a  solid  basis 
in  PBL  theory,  at  least  provides  a  departure  point 
for  interpreting  the  extremely  complex  changes  in 
momentum,  heat,  and  salt  fluxes  that  occur  across 
the  marginal  ice  zone. 
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